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1 Introduction 



Deep Inelastic Scattering (DIS) [l], fj] of electrons on protons provides the clas- 
sical test of the Quantum Chromodynamics (QCD)j3|, 0] at short distances. 
A complete knowledge of the nucleon structure (parton picture of hadrons) , 
including the gluon density, is fundamental for the design of future proton 
collider experiments, but more immediate interest in the gluon density cen- 
tres on the main role it plays for physics at low x and moderate Q 2 . With 
the advent of HERA, we are working in a new kinematical region both in Q 2 
and x. Whereas at large Q 2 the logarithmic scaling violations predicted by 
QCD are expected, the direction towards small x leads to a new kinematical 
region, where new dynamical effects are expected to occur at a significant 
level. 

The standard perturbative QCD framework predicts that at a regime of 
low values of the Bjorken variable x (x ~ 10~ 4 ) and large values of Q 2 , a 
nucleon consists predominantly of the sea quarks and gluons. The gluons 
couple only through the strong interaction, consequently the gluons are not 
directly probed in the DIS, only contributing indirectly via the g — > qq 
transition. Therefore, its distribution is not so well determined as those of 
the quarks. 

Recently, a direct determination of the gluon density was performed in 
the kinematical region x > 10~ 3 [|J. A leading order determination of the 
gluon density in the proton was performed by measuring multi-jet events 
from boson-gluon fusion in deep-inelastic scattering with the detector HI at 
the electron-proton collider HERA. The data agree with the dramatic rise 
of gluon density with decreasing x predicted by the traditional analysis of 
F 2 (x, Q 2 ) ||, ||. These analysis have performed next-to- leading order (NLO) 
QCD fits based on the Altarelli-Parisi (DGLAP)|7| evolution equations of 
the parton densities. The initial densities should be determined by the ex- 
periment at some scale Q 2 = Q^. The basic procedure is to parameterize 
the x dependence at some low Ql, but where perturbative QCD should be 
still applicable, and then evolve up in Q 2 using NLO DGLAP equations to 
obtain the parton densities at all values of x and Q 2 of the data. The input 
parameters are then determined by a global fit to the data. 

The parton densities are well determined for x > 0.5, where there are 
many different types of high-precision constraints on the individual parton 
distributions. The exception is the gluon density which is mainly constrained 
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by (i) the momentum sum rule, (ii) prompt photon production and (iii) the 
scaling violations of F 2 . On the other hand at small x (x < 10~ 3 ) we have 
only the measurement of the structure function F 2 (x, Q 2 ), and many differ- 
ent partonic descriptions at small x, for example, DGLAP[|7| and BFKL[|]. 
Therefore it is important to obtain more information in this region, mainly 
of the gluon density, which is the dominant parton at small x. 

Using the fact that quark densities can be neglected and that the non- 
singlet contribution F 2 NS can be ignored safely at low x in the DGLAP equa- 
tions, the equation for F 2 becomes, for four flavours 

dF 2 (x,Q 2 ) 10a s f 1 x fx 2 \ 

where g(x, Q 2 ) is the gluon density, xg(x, Q 2 ) = G(x, Q 2 ) is the gluon density 
momentum, a s = a s (Q 2 ) is the strong coupling constant and the splitting 
function P qg (x') gives the probability of finding a quark inside a gluon with 
momentum x' of the gluon. 
In LO (leading order) 

P qg {x') = \{x* + {l-x'Y\ . (2) 
Substituting x' = 1 — z we can write 

dF 2 ( x ,Qi)_m a , f- dzPqg{l _ z)G t * #\ _ (3) 



dlogQ 2 9n Jo 



Moreover, in LO P qg (z) = P qg (l — z), then 



dF 2 (x,Q 2 ) 10a s 



dlogQ 2 9vr 



l-x 



X 



dzP qg (z)G -,Q ) ■ (4) 



This equation describes the evolution of the proton structure function 
F 2 (x, Q 2 ) in LO at small x. From experiments one can measure the structure 
function and then determine its derivative dF ^Q2 ■ Some methods [^, |10 



were proposed in the literature to isolate the gluon density by its expansion 
inside expression (f|), but have different results for G(x, Q 2 ). We demonstrate 
here that one of the methods is not totally correct and that, when it is cor- 
rected, the results of both methods can be obtained by using different points 
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of expansion from one general expression. This article will be organized as 
follows. In sect. 2 a revision of the methods proposed in the literature is 
made. We show that one of these methods is not totally correct and we 
present a modification. In sect. 3 we propose a general expression for the ap- 
proximate determination of the gluon density and determine the appropriate 
points of expansion using data for dF ^ g q2 , obtained by Hlfnj, Zeus[[L^] and 
also for a screnning model jn| recently appeared in the literature. Finally, 
in the sect. 4, we present our conclusions. 



2 Analysis of Approximative Methods 

The methods of approximate determination of the gluon density are based 
on the simplification of the convolution P qg ®g (eq. (JJ)) by the expansion 
of the gluon density. The result is the gluon density G(k.x) proportional to 
the derivative of F 2 (x, Q 2 ), where the constant k is associated to the choice 
of a point of expansion and x = xb- In DIS we do not probe the gluon 
directly, but instead via the process g — > qq. The longitudinal fraction x g 
of the proton's momentum that is carried by the gluon is therefore sampled 
over an interval bounded below by the Bjorken xb variable xb = , where 
as usual p and q are the four momenta of the incoming proton and virtual 
photon respectively and Q 2 = — q 2 . Therefore, the choice of the point of 
expansion, and consequently of k, is associated with the relation between 
x g and xb- We discuss the methods of approximate determination of gluon 
density proposed in the literature below. 



2.1 Prytz's Approach 

Substituting P qg (x') in the expression (|J) and then expanding the gluon dis- 
tribution around z — ^, one gets || 



dF 2 (x,Q 2 ) ^ te. G{2x) [ 1 -* dz [(l-z)'+*>] + 



dlogQ 2 9tv Jo 

Approximating the upper integration limit 1 — x ~ 1 we can write 
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dF 2 (x,Q 2 
dlogQ 2 



' a °-G(2x) I' dz[(l - zf + z 2 } 
Jo 



9tt 
10a, 



G(2x) 



(6) 
(7) 



Therefore the gluon distribution can be expressed by 



G{2x) 



27tt dF 2 (x,Q 2 
10a s dlogQ 2 



(8) 



The author of ref. |9j claims that this relation is valid within 20% of accu- 
racy, by taking into account many uncertainties: experimental measurement 
of the Q 2 -dependence of F 2 , the value of a s and the discarding of any typi- 
cal low-x effects calculated beyond the DGLAP approximation, for example 
recombination effects [P^]. Estimates of gluon density M have demonstrated 
that this approximation agrees with the gluon density obtained with other 
indirect determinations made by HI and Zeus Collaborations. 



2.2 Bora's et al. Approach 

A different method is based in expansion of gluon distribution around z 



Substituting the result in the expression (H) one gets [10 



dF 2 (x,Q 2 ) 5a, 
dlogQ 2 



9tt 



9tt 



G{x) 



2(1 - x) 



(1-x) 2 + (l-x) + 



+ 5a s f(l-x) 2 2(1 -xf (1 -x) 4 




(9) 



Therefore 




where 



dF 2 (x,Q 2 ) ^5a SA ^i G{x) + m^ x dG) (lQ) 



A( x ) = 2(1 x)3 - (1 - xf + {l-x) , (11) 
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and 



B(x) = {l ~ X? - 2{l ~ xf + (1 ~ g)4 



(12) 



Bora et al. [10 claim that this expression can be approximated by 



dF 2 (x,Q 2 ) 5a s [A(x) + B(x)f 



G x 



B(x) 



dlogQ 2 9tt A(x) + 2B(x) \ A(x) + B(x 
B(x) \ 9tt A(x) + 2B(x) dF 2 (x,Q 2 



x 



G \ x + 



-x 



v A(x) + B(x)~J ' 5a s \A(x) + B{x)} 2 dlogQ 2 
In the limit x —>■ the last equation becomes 



(13) 
(14) 



r:( '4 \ 9n4 dF 2 (x,Q 2 ) 
,3 / 5a s 3 dlogQ 2 



(15) 



In figure (P we compare the results of dF ^ g Qi - for the approximate re- 
lations of Prytz (0), Bora et a/.(|T3"D and the complete expression (f|) using 
the gluon distribution GRV94(LO) M. In our calculations the value of a s is 
obtained at = 0.232GW 2 . The fi gure shows that Prytz approximation 
is more accurate. 

The result ([TJJ) differs from Prytz result (Q). Bora et al. explain that the 
difference arises because they have retained the x dependence in the upper 
limit of integration in equation (f|) after the expansion of the gluon density. 
In order to show this is not the case, we have computed the expression (||) 
without approximation in the upper integration limit. One gets 



dF 2 (x,Q 2 
dlogQ 2 



5as 
~9tt 



[(l_ a; )_(l_ a ;)2 + _( 1 _ a .)3 

[l-xf -\{l-xf + \(l-x) 4 



1-x) 



x) 2 + j(l 



X) 



G{2x) + 
1 

" 2 



dG(z 



dz 



(16) 



Consequently, 
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<iF 2 (x,Q 2 ) 5a s . .„ , , 9 2._ ,o. 



i(l-x)-|(l-x)_ 2 + i(l-.r 

2 
3 

In the limit x — > 0, this equation becomes 



G * + 2 i-ri- l)+ a-,p • (17) 



(19) 

This result is identical to the Prytz one, demonstrating that Bora's argu- 
ment is not correct. 

The Bora's method requires that 



A(x)[A(x) + 2B(x)} 
A(x) + B(x) 
A(x) + 2B(x) 



, (20) 

1 . (21) 



In figure (0) we present the behaviour of these expressions at small x. 
We immediately see that for small x these approximations are not valid. 
Consequently the Bora approximation is not valid in this x range. 

The more adequate approximation of expression fllCf) is 

dlogQ 2 9tt V ; V M x ) J 

In the limit x — > , the equation becomes 



<«F,(x,Q») 6*2 /3. 



d/o#Q 2 9tt 3 V2 
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In figure ([3]) we presented the behaviour of Prytz (|7|) and Bora ( |T3| ) results 
and one obtained from expression (|23|) , called Bora corrected. The expression 
(53) also differs from the expression of Prytz (|7p. This is due to the election 
of a different point of expansion of the gluon density in both cases. In the 
next section we demonstrate that both results can be obtained from a general 
equation by choosing the adequate point of expansion. 



3 Expansion at an arbitrary point of expan- 
sion 

Using the expansion of the gluon distribution G (jzj) at an arbitrary z = a 
and retaining terms only up to the first derivative in the expansion, we get 



dlogQ 2 ~ 9tt 1 3^ 



.G 



x 



a + 



f(l-*) 3 -f(l 



;i -xy + (i-x 

x) 2 + ±(l-x)X 



l-oc\ f (1 - x) 2 - (1 -x) + 1 

In the limit x — > , the equation becomes 



(24) 



dF 2 (x,Q 2 
dlogQ 2 



5a s 2 



-G 



9tt 3 



x 



3 

a 

a \2 



When the points a= \ and a = are used, we get 



(25) 



dF 2 (x,Q 2 ) 
dlogQ 2 

dF 2 (x,Q 2 ) 
dlogQ 2 



5 -^ 2 -G(2x) 
9tt 3 V ' 

5a,2_/3 



(26) 
(27) 



Therefore the equation ( p^) can be reduced to Prytz (^) and Bora corrected 
(|23|) , in their respective points of expansion. 

The gluon distribution for a point of expansion a < 1 can be expressed 

by 
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G 



x /3 



I- a \2 



a 



9tt 3dF 2 {x,Q 



2\ 



5a s 2 dlogQ 2 



(28) 



In figure we present the results of dF ^Q 2 - at some points of expansion 
a using the gluon distribution GRV94(LO), which fits better the recent data 
of HI Collaboration fllS] . At the current data of dF ^ g Q2 the points a > 0.5 



are favoured. This means that in this kinematical region the longitudinal 
momentum of the gluon x g is more than twice the value of the longitudinal 
momentum of the probed quark (or antiquark) in DIS. 

One can ask if the approximative determination of the gluon density from 
F 2 (x, Q 2 ) scaling violations can indicate the presence of new effects at low-x, 
for example, recombination. In figure (|5[) we present the behaviour obtained 
using fl2"E| ) for distinct points of expansion and dF ^Q2 ^ data calculated using 



the Glauber(Mueller) approach ||13||. This approach includes shadowing cor- 
rections and dF ^ g q2 ^ can be calculated using the total cross section of the 
gluon pair with the nucleon in the eikonal approach. 

We can conclude that, when compared with the behaviour of GRV94(LO), 
the screnning can be cancealed for some values of a and that the better 
choices (with greater sensitivity to screnning) are in the range 0.5 < a < 0.8. 
This conclusion agrees with that of Ryskin et al. |18j that estimate the value 
of the longitudinal gluon momenta x g in approximately 3 times larger than 
the Bjorken xb, that corresponds to the expansion at a = 0.75. 



4 Conclusions 

The gluon is by far the dominant parton in the small x regime. Its distri- 
bution is not well determined as those of the quark, and there are several 
methods proposed in the literature in order to determine it in different re- 



gions of x. More recently, new methods 0, [10], [16], [17]] also based on the 
dominant behaviour of the gluon distribution at small- a; were proposed. 

In this work we discuss the approximative determination of the gluon 
density at \ow-x from F 2 (x,Q 2 ) scaling violations 0, [TtJ. We demonstrate 
that one of these methods is not correct. Moreover we proposed one general 
expression for approximative determination of the gluon density at arbitrary 
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point of expansion of G [jz^j ■ Using data of HI and Zeus Collaborations 
n, [12] and results obtained from Glauber (Mueller) approach |13[] that in- 
cludes shadowing corrections, we conclude that the more suitable points of ex- 
pansion (with greater sensitivity to screnning) are in the range 0.5 < a < 0.8. 
This conclusion agrees with recent results obtained by Ryskin et al. [18| . 

The approximative determination of the gluon density obtained in this 
work was developed in the DGLAP framework. Consequently the approxi- 
mation is valid in the DGLAP limit, which surely breaks down at sufficiently 
small x. When aJn- 1 we should resum also the aJn- contributions. In 

a X a X 

LO, (a s ln^) n resummation is acomplished by the BFKL equation ||. There- 
fore, a rigorous determination of the transition region to small-x dynamics 
is very important in order to isolate correctly the gluon distribution. 
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Figure Captions 



Fig. [T]: Results of the ^gg 2) obtained from relations of Prytz (0), Bora (0) 
and (|), using GRV94(LO)§. Data from HI JTTJI and Zeus[|lJ at Q 2 = 20 
GeV 2 . 



Fig. 0: Behaviour of approximations used by Bora et al. [FLO . 



Fig. ^ Same as for figure ([I]), including results of equation ( |2"BD presented 
for Q 2 = 20 GeV 2 . 



Fig. [|: Behaviour of dF ^ g q2 ~ at several points of expansion a. Also 
plotted the behaviour obtained using GRV94(LO) and the expression ([|) at 
Q 2 = 20 GeV 2 . 



Fig. U Gluon distributions obtained of ( p8|) using dF ^ g Q2 ^ with screnningp3 
Also shown the gluon distribution of GRV94(LO). 
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Figure 3: 
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Figure 4: 
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